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INTRODUCTION 


The  year  1980  marks  the  centenary  of  the  discovery  of  piezoelectricity 
by  the  brothers  Curie,  and  the  fifty-ninth  year  of  the  quartz  oscillator. '* 
The  past  thirty  years  or  so  have  seen  the  development  of  cultured  quartz 
from  coin-sized  pieces,  shown  in  Fig.  1,  to  large,  multi-kilogram  crystals 
that  constitute  the  bulk  of  the  usage  today  (Fig.  2).  Since  the  introduc¬ 
tion  of  the  quartz  crystal -control led  oscillator  by  Cady,  each  decade  of 
time  has  seen  approximately  one  order  of  magnitude  frequency  stability 
improvement.  Surrounding  the  development  of  high-stability  oscillators  in¬ 
corporating  crystal  vibrator  control  has  been  the  on-going  interplay  between 
improvements  in  resonator  technology  and  causes  of  ever-smaller  perturba¬ 
tions  in  resonator  frequency  and  equivalent  circuit  parameters. U42* Progress 
in  timekeeping  is  shown  on  a  broad  time  scale  in  Fig.  3,  and  a  comparison 
of  crystal  resonator  performance  with  atomic  and  molecular  frequency  stan¬ 
dards  is  shown  in  Fig.  4;  both  Figures  are  adapted  from  Reference  2.* 

Further  improvements  will  depend  upon  less-obvious  factors  than  have  been 
considered  to  date.  This  technical  report  treats  a  number  of  such  topics. 

MAGNETIC  FIELDS 

Quartz,  and  other  high  purity  dielectric  crystals,  can  be  considered 
to  be  nonmagnetic  for  our  purposes  here.  A  vibrating  crystal  with  metallic, 
nonmagnetic  electrodes  is,  however,  affected  by  external  magnetic  fields. 
Eddy  currents  are  produced  in  the  electrodes  by  their  motions  in  the  fields; 
these  change  the  critical  frequencies  of  the  vibrator,  and  lower  its  Q.  The 
situation  considered  is  that  of  a  resonatot  vibrating  in  a  pure  shear  mode 
with  applied  static  magnetic  field  parallel  to  plate  thickness  as  depicted 
in  Fig.  5.  The  resonator  leads  are  attached  to  an  active  network  comprising 
the  remainder  of  the  oscillator,  and  the  combination  adjusted  in  frequency 
so  that  the  vibrator  operates  at  zero  reactance  at  all  times. 

For  a  pure  shear  mode,  the  retarding  force  density  at  the  plate  sur¬ 
faces  is 


-  j  u  o  B2  u,  (1) 

with  a  the  electrical  conductivity  of  the  electrode,  and  u  the  mechanical 
displacement  at  the  surfaces.  Electrodes  of  thickness  h'  will  produce 
resistances  of 


R  =  A  B2  h’  (2) 

3* 

for  areas  of  value  A.  Using  conventional  notation  ,  this  leads  to  an  ex¬ 
pression  for  the  imput  admittance  of  the  single  mode  vibrator,  including  the 
effects  of  electrode  inertia4  : 


*See  list  of  references  beginning  on  page  72. 
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OF  OSCILLATORS  AS  FUNCTION  OF  SAMPLE  TIME. 


Yin  =  j  w  C  +  j  uj  Cf). 

(3) 

{ (X/(k‘  tan  X)  -  1  -  X'  u/k  )  +  ( jXRY  ,  /k  )r 

Resonance  and  anti  resonance  frequencies  are  determined  from 
the  susceptance  function  as  the  roots  of 

the  zeros  of 

(X-(k‘  +  u  X  )  tan  X)  •  (1-  u  X  tan  X) 

+  X  '  (R  Y0)2  •  tan2  X  =  0. 

(4) 

When  R=0,  the  first  expression  in  brackets  yields  the  condition  for 
resonance,  and  the  second  that  for  antiresonance.  From  the  expression  for 
Yin  The  network  in  Fig,  6  is  obtained,  with  the  electrode  inertia  represent¬ 
ed  by  the  inductances 

Li  =  L2  =  L  =  p  h  u  A. 

(5) 

The  network  components  are  labeled  with  their  physically  corresponding  ef¬ 
fects.  When  R  is  not  zero,  the  inductances  are  replaced  by  series  R-L 
combinations,  and  when  electrode  mass-loading  is  neglected,  the  circuit 
appears  as  given  in  Fig.  7.  Figure  8  is  the  bisected  version  of  Fig.  7  with 
massless  electrodes  (u=0). 

Departures  from  the  resonance  frequency  as  function  of 
oerturbation,  with  RY0  the  small  parameter.  The  expression 
shift  is 

R  are  found  by 
for  frequency 

Af/f  »  -  ( 2/ tt ) 2  (RY0)2, 

(6) 

with  RYo  =  o  h'  BZ  /  (dv). 

(7) 

i  the  quartz  mass  density,  and  v  the  acoustic  velocity  of  the  mode  in  ques¬ 
tion.  It  is  seen  from  (6)  and  (7)  that  the  frequency  shift  is  quartic  in 
the  magnetic  field  strength.  Figure  9  gives  frequency  shift  values  for  an 
AT  cut  having  10,000  A  copper  electrodes:  only  for  fields  about  ten  thousand 
times  larger  than  that  of  the  earth's  (the  first  entry)  is  the  effect  per¬ 
ceptible. 

This  analysis  has  been  extended  to  the  multimode  case,  and  to  the  sit¬ 
uation  of  a  lateral  magnetic  field,  with  corresponding  network  representa¬ 
tions. 

RESONATOR  ASYMMETRIES 

In  treatments  of  resonator  electroding,  it  is  usually  tacitly  assumed 
that  both  electrodes  are  of  equal  size  and  thickness.  The  practical  real¬ 
ity  is  usually  just  the  opposite.  It  is  very  difficult  to  guarantee  that 
the  thicknesses  are  equal  in  the  first  place;  to  make  matters  worse,  final 
frequency  adustment  is  usually  made  with  a  burst  of  material  applied  to  a 
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FIGURE  6.  AnalO'i  Network  of  Resonator  with  Electrode  Imbalance 


CRYSTAL  PLATE 


FIGURE  7.  Analog  Network  of  Resonator  in  Magnetic  Field 


single  side.  Add  to  this  situation  the  facts  that  registration  of  both 
electrodes  is  never  perfect,  that  the  electrode  tabs  are  asymmetrically 
disposed,  and  that  sometimes  a  larger  electrode  is  used  on  one  side  to  im¬ 
prove  statistical  unii  rrnity  of  the  motional  capacitance  by  compensating 
for  misregistration,  t 'en  it  is  seen  that  symmetry  of  the  plate  resonator 
about  the  mid-plane  is  seldom,  if  ever,  achieved. 

The  single-mode  case  is  treated,  by  way  of  example,  for  a  simple  elec 
trode  mass  imbalance  in  che  one-dimensional  approximation.  The  equivalent 
network  is  that  of  Fig.  6  with  the  mechanical  transmission  line  ports 
attached  to  unequal  transmission  lines  representing  the  thick  electrodes. 
When  the  electrodes  are  thin  enough  for  propagation  of  the  acoustic  waves 
within  the  electrodes  to  be  neglected,  the  transmission  lines  are  replaced 
with  lumped  inductances;  the  mass  imbalance  results  in  unequal  inductance 
values  Li  and  L2,  as  shown,  and  prevents  the  network  from  being  bisected, 
as  in  Fig.  8.  When  bisection  can  take  place,  the  transmission  line  repre¬ 
senting  the  crystal  plate  of  thickness  2h  is  replaced  by  one  of  half  this 
value  and  the  construction  for  the  critical  frequencies  is  that  given  in 
Fig.  10;  the  formulas  are  those  of  the  bracketed  terms  in  (4)  separately 
set  to  zero.  Figure  11  gives  the  normalized  frequency  variable 

-  (■/!')  (fRiM)/fA0U> 

versus  piezoelectric  coupling  k  for  the  first  three  harmonics  M=l,  3, and  5 
An  AT-cut  resonator  with  symmetrical  electrodes  of  equal  thickness  pro¬ 
duced  the  mode  spectrograph  at  the  top  of  Fig.  12.  Strong  responses  are 
found  at  the  odd  harmonifs,  and  nothing  is  visible  at  the  even  harmonics. 

It  is  appropriate  now  to  quote  from  Cady1,  p.  308: 

"There  are  circumstances  under  which  a  plate  may  conceivably 
vibrate  in  a  thickness  mode  at  or  close  to  an  even  harmonic 
frequency.  The  shape,  size  and  location  of  the  electrodes 
may  be  such  as  to  produce  a  driving  field  in  the  plate  that 
varies  in  the  direction  of  the  thickness.  Or  the  plate  may 
be  twinned  or  have  other  defects  such  that  the  excitation 
is  not  uniform.  Finally,  the  plate  may  be  in  contact  with 
an  electrode  of  considerable  mass,  so  that  in  effect  one 
has  a  composite  resonator." 

When  the  electrodes  are  of  unequal  mass,  and  circuit  bisection  is  not 
possible,  the  construction  becomes  that  of  Fig.  13,  where  the  critical 
frequencies  are  determined  by  intersection  with  a  tangent  function  of 
argument  2X  instead  of  X.  The  expression  for  input  impedance  is 

k?  (tan  X/X) 

Zin  =  (1/jX,,)  •  U - 1  ,  (8) 

a2X2 

((l-HX  tanX)  + - •) 

( 1+pX  cotX ) 
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FIGURE  10.  Construction  for  Resonance  and  Antiresonance  Frequencies 


FIGURE  11.  Normalized  Resonance  Frequency 


FIGURE  13.  Construction  for  Anti  resonant  Frequencies  with  Imbalance 


where  '  =  +  u,)/2  and  a  =  (n  -  uJ/2.  Eqtation  (8)  has  been  written 

in  terms  of  argument  X  for  the  tangent  and  cot.  ngent  functions.  The  <inti- 
resonant  frequencies  (denoted  f.(M)  for  M  odd,  and  f>(M)  for  M  even)  are 
found  from  the  vanishing  of  theAbracketed  term  in  (Hf‘  as  roots  of 

-(  M  t  *  v1  -  )  X 

tan  2X  =  '  ~  •  19) 

(1  -  n,  a,  X‘  ) 

These  are  shown  in  Fig.  13  along  with  the  high  and  low  frequency  asymptotes 
of  the  right  hand  side  of  (9).  The  roots  for  the  resonant  frequencies  are 
found  from  the  zeros  of  (8). 

The  even-harmonic  resonances  produced  by  the  unequal  mass-loadings 
on  the  crystal  plate  surfaces  (pj  t  u,)  are  very  sharp  compared  to  the  ode 
harmonics  as  may  be  seen  from  Fig.  14,  which  is  drawn  for  a  case  with  large 
piezo-coupling  and  large  mass  asymmetry  in  order  to  widen  the  resonance 
regions  of  both  the  even  and  odd  harmonics.  It  is  thus  possible  to  use 
such  resonances  for  stabilizing  the  frequency  of  an  oscillator  -  at  least 
in  principle.  Figure  15  shows  a  blow-up  of  the  second  harmonic  of  the  case 
of  Fig.  14;  Fig.  16  shows  the  fourth  harmonic,  and  Fig.  17  shows  the  four¬ 
teenth  harmonic.  Compare  the  width  of  the  resonance  in  Fig.  17  with  that 
of  Fig.  18  showing  the  fifteenth  harmonic.  In  general,  for  large  M,  the 
normalized  frequency  of  response  approaches  (Ml). 

When  the  resonator  shown  in  the  top  of  Fig.  12  was  overcoated  with 
electrode  material  on  one  side  only,  the  mode  pectrum  shown  in  the  bottom 
of  the  figure  was  produced.  The  second  harmon  c  becomes  present  due  to 
the  imbalance,  and  it  is  sharp;  however,  it  is  also  very  weak.  Figure  19 
shows  a  photograph  of  the  spectrograph  given  i  i  the  top  of  Fig.  12;  Fig.  20 
is  a  photograph  of  the  spectrograph  for  the  case  of  mass  imbalance  to  the 
same  scale  as  Fig.  20;  Fig.  21  shows  a  magnified  view  of  the  photograph  of 
Fig.  20;  it  corresponds  to  the  lower  drawing  in  Fig.  12. 

Another  unsymmetrical  structure  that  lead1,  to  even  harmonics  is  the 
composite  resonator;  (see  Reference  5  for  additional  references).  This 
structure  is  shown  in  Fig.  22,  with  the  equivalent  network  given  in  Fig. 

23.  Here  the  quartz  substrate  must  be  represented  by  a  transmission  line 
because  of  its  thickness;  as  it  shrinks  in  size  the  situation  of  Fig.  6 
with  inductors  is  approached. 

Apart  from  producing  even  harmonics,  the  presence  of  asymmetry  in 
resonators  (due  to  unequal  electrodes,  bevelling,  contouring,  etc.)  is 
very  probably  a  major  cause  of  "activity  dips,'  or  "bandbreaks . "6"8  These 
are  illustrated  in  Fig.  24,  which  is  taken  from  Reference  8.  Both  the  fre¬ 
quency  versus  temperature  and  the  conductance  versus  temperature  curves 
suffer  anomalous  dips  as  the  interfering  mode's  influence  is  manifest  at  a 
particular  temperature.  The  second  harmonic  mode  "divided  by  two"  via  the 
nonlinear  elastic  constants9  can  produce  dips  at  the  fundamental. 
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FIGURE  16.  Spectrograph  Near  the  Fourth  Harmonic 


FIGURE  17.  Spectrograph  Near  the  Fourteenth  riarir'onic 
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FIGURE  19.  Photograph  of  Spectrum  of  Balanced  Resonator 
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FIGURE  21.  Magnification  of  Spectrum  of  Figure  20 
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COMPOSITE  RESONATOR  STRUCTURE 


L 


UPPER 

ELECTRODE 


THIN-FILM 

DRIVER 

LAYER 


LOWER 

ELECTRODE 


QUARTZ 

SUBSTRATE 


FIGURE  23.  EQUIVALENT  NETWORK 

OF  COMPOSITE  RESONATOR 


TEMPERATURE  (  °C  ) 


ELECTRODE  INFLUENCES 


It  is  well-known  that  the  frequency-temperature  (f-T)  curve  of  an  AT 
plate  resonator  measured  sufficiently  slowly  assumes  a  cubic  shape,  as 
depicted  in  Fig.  25.  Small  changes  in  orientation  angle  •>,  for  rotation 
about  the  crystallographic  Xi  axis,  produce  slope  changes  at  the  inflec¬ 
tion  point  of  -5.08  x  10"(7K,  °  o.  A  variety  of  other  influences  lead  to 
nonnegl igible  slope  variations  as  well.  Figure  25  shows  the  effect  of 
placing  a  load  capacitor  in  series  with  the  vibrator  to  produce  frequency 
fL,  as  is  often  done  for  frequency  trimming.  The  slope  shift  reduces  the 
peak-to-peak  frequency  deviation  and  must  be  taken  into  account  in  design 
of  temperature  compensated  crystal  oscillators  (TCXOs);  it  comes  about 
because  of  the  temperature  coefficient  of  the  piezoelectric  coupling  fac¬ 
tor.10  A  similar  effect  occurs  when  a  resonator  is  operated  at  different 
harmonics,  and  for  the  same  reason. 

The  mass  of,  and  intrinsic  stress  in  electrode  coatings  also  affects 
the  f-T  curve.  Some  time  ago  the  effect  of  equal  electrode  coatings  was 
demonstrated  in  regard  to  frequency  shift  and  rotation  of  the  static  fre¬ 
quency-temperature  curve  of  a  resonator.11-13  The  calculations  dealt  only 
with  the  influence  of  inertial  mass;  no  elastic  properties  were  assumed 
for  the  electrodes,  nor  were  they  presumed  to  apply  any  stresses  to  the 
crystal.  This  was  done  so  that  the  inertial  influence  of  the  electrode 
could  be  subtracted  from  elastic  and  stress  effects.  Results  from  the  one 
dimensional  theory  for  the  change  in  apparent  orientation  angle  are  given 
in  Fig.  26  and  Fig.  27  as  function  of  harmonic  number  M  and  mass  loading  . 

A  comparison  of  apparent  orientation  angle  shifts  as  function  of  mass 
loading  and  harmonic  number  is  given  in  Fig.  21  for  AT  plates  of  5MHz  fun¬ 
damental  frequency,  and  plate/electrode  diameters  of  14  and  6  mm,  respec¬ 
tively,  coated  with  evaporated  aluminum,  sputtered  chromium,  and  sputtered 
indium.  The  electrodes  are  of  equal  thicknesses.  It  is  seen  that  A1  and 
In  have  slopes  that  agree  with  inertial  mass  theory,  but  are  displaced;  the 
curves  for  Cr  are  consistent  among  themselves,  but  differ  in  sign  from  the 
theory.  The  results  in  Fig.  28  indicate  the  importance  of  electrode  ma¬ 
terial,  deposition  method,  and  conditions  of  deposition  on  the  behavior  of 
the  finished  resonator.  Figure  29  gives  the  experimental  data  for  alumin¬ 
um  evaporated  electrodes  as  function  of  p  and  M.  Figure  30  shows  the 
decrease  in  inflection  temperature  with  increasing  mass  of  evaporated  alu¬ 
minum,  for  resonators  with  the  parameters  described  above. 

It  should  be  borne  in  mind  that  stresses  in  the  electrodes  can  also 
produce  activity  dips,  even  when  the  electrodes  are  of  equal  mass.  Conse¬ 
quently,  the  practice  of  furnishing  AT  quartz  plates  with  a  thin  layer  of 
Ni  or  Cr  to  insure  adhesion  of  gold,  e.g.,  ought  to  be  re-examined  in  view 
of  the  stress  levels  produced  hy  these  metals,  both  intrinsically,  and  as 
a  function  of  temperature. 
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figure  29.  APPARENT  ANGLE  SHIFT  IN  MINUTES  OF  ARC 
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ANGULAR  MISORIENTATION 

Subtle  influences  on  the  static  frequency-temperature  curve  were 
briefly  discussed  in  the  previous  section  in  terms  of  apparent  orientation 
angle  changes  in  0,  for  AT  resonators.  The  most  general  class  of  thickness 
mode  plate  resonators  is  the  family  of  doubly  rotated  cuts.14  These  are 
compared  in  Fig.  31  with  singly  rotated  cuts,  of  which  the  AT  cut  is  the 
most  often  used.  Also  shown  in  the  figure  are  the  loci  of  zero  temperature 
coefficient,  and  the  angular  coordinates  of  popular  quartz  cuts.  The  line 
$  =  0  is  a  line  of  symmetry,  so  that  singly  rotated  cuts  are  only  quadrat- 
jcally  sensitive  to  misorientations  in  <j> .  For  the  AT  cut,  misorientations 
-A 4  may  be  compensated  for  by  changes  A6,  obtained  from 

AO  =  -1.77  x  lO-3  (a*)2, 

with  A<t>  in  minutes  of  arc,  and  Ae  in  seconds  of  arc;  this  relationship  is 
portrayed  in  Fig.  32. 

X-RAY  ANGLE  CHANGES1'1 

The  angles  locating  the  X-ray  planes  in  a  crystal  are  functions  of 
temperature  and  external  stress.  For  example,  the  plane  at  38°  13'  used 
for  orienting  AT  cuts  is  only  found  at  this  angle  at  a  certain  temperature. 
For  high  precision  applications  the  angle  shift  with  temperature  is  not 
negligible.  We  give  an  approximate  calculation  to  establish  the  order  of 
magnitude  of  the  angle  change. 

Consider  a  rotated-Y-cut  plate  of  crystal  in  the  rhombic,  tetragonal, 
trigonal,  hexagonal,  or  cubic  systems.  (These  have  no  off-diagonal  thermo¬ 
elastic  coefficients.)  At  temperature  T  the  plate  angle  will  be  deternrnpd 
from 


tan  e  =  (X2/X3).  (10) 

With  a  change  in  temperature  to  T+aT,  the  axes  become,  in  the  linear 
approximation 


X2  =  X2  (1  +  a22  AT) 

Xi  =  X  3  (1  +  CX33  AT) 

We  are  thus  led  to  a  relation  between  Ae  and  aT: 

tan  Ae  =  { ( A- 1 )  tan  0  /  (1  +  A  tan2  0  )} 


(11) 


(12) 


where 

A  =  (1  +  a22  AT)/(1  +  a33  AT).  (13) 
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Doubly  Rotated  Cuts  and  Zero  Tem^ra 


±  A<X> 


FIGURE  32.  Delta  Theta  versus  Delta  Phi  tor  AT  Cuts 


<■»  /• 


As  a  consequence  of  (12)  and  (13)  we  have 

A6  =  '2  (au  -  a33)*  AT  •  sin  28.  (14) 

Figure  33  gives  the  X-ray  angle  change  with  temperature;  beyond 
about  10  kelvins  the  linear  approximation  is  not  accurate  except  as  to 
order  of  magnitude.  For  the  tetragonal,  trigonal  and  hexagonal  systems, 
since  a?2  =  au,  the  result  is  independent  of  angle  <p.  In  particular,  for 
the  upper  zero  temperature  locus  for  thickness  mode  plates  f,  is  approxi¬ 
mately  constant  and  the  result  holds  therefore  for  the  AT,  FC,  IT,  and  SC 
cuts  (cf.  Fig.  31).  Relation  (14)  shows  that  ap  depends  on  the  difference 
between  an  and  u33  for  rotated-Y-cuts  (or  doubly  rotated  cuts  when  n:,y  = 
an).  For  the  cubic  system  this  difference  is  zero,  so  ao  is  always  zero 
in  this  system.  For  quartz,  the  effective  value  of  a  in  the  thickness 
direction  of  rotated-Y-cut  plates  is  shown  in  Fig.  34;  it  is  given  by 

a  =  1U  COS2  6  +  a  sin2  8;  (15) 

for  the  Y  cut,  a  =a0o  =  a, ,=  13.72  x  10-6/K,  while  for  the  Z  cut  (6  =  -90), 
a  =a33  =  7.48  x  10'6/K. 

TWINNING 


The  types  of  twinning  that  may  occur  in  quartz  are  shown  in  Fig.  35 
along  with  the  changes  in  axes  and  handedness.16’1  ’  We  briefly  explore  the 
effects  of  twinning  on  the  piezoelectric  constants  of  crystals  in  Class  18, 
which  includes  quartz.  For  the  three  types  of  twinning  considered,  the 
operations  involved  are 

(a)  Dauphine:  Xj-^-Xj;  X2-*~X2;  X  -*-+X3 . 


(b)  Brazi  1 : 


en-e. 


2  2 


e  -»+e 

14  C1 


X1-v-X1 ;  X2^+X2;  X  ++X, 


ell!+_eii;  el4^ei<.  • 

(c)  Combined:  X^+Xj;  Xo-*—X2 ;  X  ++X3. 

eir+eni  e14^ei4- 

Consider  e'26,  the  piezoelectric  constant  driving  the  pure  shear  mode 
in  rotated-Y-cuts  of  quartz: 


Untwinned:  e'26  =  (en  cos  e  +  e14  sino)  *  (-cose). 

(a)  Dauphine:  e'26  =  (-en  cose  +  e;4  sine)  •  (-cose). 

(b)  Brazil:  e'26  =  (  en  cose  +  e14  sine)  ■  (+cose). 

(c)  Combined:  e'26  =  (  eu  cose  -  e|4  sine)  ’  (-cose). 
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FIGURE  35.  El 


Consider  further  a  plate  of  unity  area,  with  a  fraction  "a"  of  this 
area  twinned.  For  simplicity,  assume  the  twinning  is  complete  throughout 
the  thickness,  and  that  the  distribution  of  motion  is  uniform  even  though 
twinned  (although  we  keep  in  mind  Cady's  remark,  above).  Then  the  compos¬ 
ite  e'26Should  be: 

(a)  Dauphine: 

e '”26  =  ( ( 1  - 2a )  en  cose  +  e  sine)  ‘  (-cose). 

(b)  Brazil : 

e'26  =  (l-2a)  en  cose  +  ( 1  -2a )  sine)  ’  (-cose). 

(c)  Combined: 

e'26  =  (en  cose  +  (l-2a)  e14  sine)  ‘  (-cose). 

In  the  case  of  Brazil  twinning: 

e"2fj  (twinned)  =  (l-2a)  '  e'2b(untwinned) .  (16) 

Equation  (16),  when  written  in  terms  of  the  resonator  capacitance  ratio 
becomes : 

r  (twinned)  =  r  (untwinned)  /  (l-2a)2.  (17) 

DYNAMIC  THERMAL  EFFECTS 

It  has  been  known  for  many  years  that  abrupt  temperature  changes 
produce  frequency  changes  in  quartz  vibrators  that  are  unpredicted  by  stat¬ 
ic  f-T  characteristics  such  as  those  of  Fig.  25.  The  dynamic  thermal  ef¬ 
fect  has  been  given  a  phenomenological  explanation  by  a  modification  of  the 
cubic  curve  to  include  time-dependent  terms  that  take  into  account  thermal 
gradients  in  the  resonator  and  support  structure.18  The  classical  f-T 
equation 

AJF  =  a0  AT  +  b0  AT  +  cQ  aT3  (18) 


is  augmented  by  the  addition  of  the  term 

(It  AT  +  $)  •  dT/dt  (19) 

A  temperature  ramp  (dT/dt  =  constant)  applied  to  the  resonator  both  ro¬ 
tates  the  f-T  curve,  through  2r,  and  translates  it,  by  means  of  a.  Figure 
36  defines  the  desired  reference  temperature  (which  may  be  at  either  zero- 
slope  point)  at  which  the  oven  ought  to  be  set,  along  with  the  actual  oven 
setting  point  achieved,  and  the  oven  cycling  range.  The  cycling  range  is 
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FIGURE  36.  Defiritions  Relating  to  Qveni/eii  Oscillators 


a  measure  of  the  quality  of  the  crystal  oven.  The  result  of  sinusoidal 
temperature  variations  of  magnitude  aT  =  60K  about  the  inflection  temper¬ 
ature  of  the  static  AT  cut  f-T  curve  is  shown  in  Fig.  37.  Temperature  var¬ 
ies  as  sin  L.>t ;  various  values  are  shown.  In  Fig.  38  is  given  the  re¬ 
sult  of  variations  for  aT  =  10K  about  the  upper  turning  point.  Figure  39 
gives  the  result  when  the  temperature  magnitude  is  reduced  to  5mK.  In  Fig. 
40  the  effects  of  aT  =  both  and  1  mK  on  the  orbit  shapes  is  seen.  Fig¬ 
ure  41  shows  the  result  of  sinusoidal  variations  in  temperature  of  magni¬ 
tude  -;a  mK  about  the  turning  point  of  the  static  f-T  curve  of  Fig.  25. 

With  a  period  of  8  hours,  the  frequency  change  is  about  3  x  10";: ;  when  the 
period  is  lengthened  to  1  week,  the  change  is  still  about  1  x  10"n,  where¬ 
as  the  static  f-T  curve  would  predict  a  change  of  only  a  few  parts  in  10"  \ 
irrespective  of  cycling  rate.  Since  thermal  transients,  temperature  cy¬ 
clings,  and  fluctuations  cannot  be  entirely  avoided,  the  dynamic  tempera¬ 
ture  effect  is  a  very  important  consideration  for  high  precision  resonators. 
The  magnitude  of  the  effect  can  fortunately  be  reduced  by  recourse  to 
doubly  rotated  quartz  cuts,  specifically  the  SC-cut.  l'*» 19 


AMBIENT  VIBRATIONS 


Insensitivity  to  acceleration  is  one  of  the  most  prominent  criteria 
by  which  to  judge  high  stability  resonators.  At  present,  the  rule-of- 
thumb  number  relating  frequency  shift  to  applied  acceleration  is  2  x  10"1' 
per  q,  although  efforts  are  underway  to  reduce  this  value.  (See  references 
given  in  Reference  20.)  It  appears  that  a  sensitivity  of  1  x  I0"n/g  is 
a  reasonable  goal  for  the  next  few  years.  These  two  sensitivity  figures 
have  been  used  in  conjunction  with  published  values  for  acceleration  levels 
in  various  environments-'1”- *'  to  produce  the  frequency  shift  entries  shown 
in  Fig.  42  and  Fiq.  43. 

From  Fig.  42  and  Fig.  43  it  is  seen  that  the  decrease  in  acceleration 
sensitivity  to  parts  in  10!1  per  g  and  beyond  is  highly  desirable.  The 
smallest  entry  in  the  figures  is  due  to  microseisms.22'26  Figure  44,  ta¬ 
ken  from  Reference  24,  shows  the  very  interesting  spectrum  of  microseisms 
observed  at  the  earth's  surface.  It  appears  that  the  spectrum  of  displace¬ 
ment  versus  period  plotted  on  equal  log-log  scales  consists  of  portions  of 
+2  and  -2  slope  peaking  rather  sharply  at  a  period  of  about  7  seconds.  The 
maximum  amplitude  of  acceleration  is  approximately  4  micro-g's,  and  is  con¬ 
stant  along  the  +2  slope  line.  For  oscillators  with  high  acceleration 
sensitivity  this  source  of  disturbance  would  have  to  be  eliminated  by  mount 
ing  supports  even  for  isolated,  underground  installations.  Building  vibra¬ 
tions  are  far  greater  than  microseisms,  and  their  noise  spectra  must  be 
considered  in  the  installation  of  precision  oscillators.  In  certain  loca¬ 
tions,  the  accelerations  due  to  earthquakes  must  be  allowed  for,  although 
the  vibrations  found  in  urban  environments  produce  accelerations  associated 
with  earthquakes  of  magnitude  3-4  on  the  Richter  scale.21  This  scale  is 
defined,  for  shallow  shocks,  as  the  1 og t n  of  the  maximum  trace  amplitude  in 
micrometers  with  which  the  standard  short-period  torsion  seismometer  would 
register  that  earthquake  at  an  epicentral  distance  of  100km. 

Two  crystal  resonator  structures  having  decreased  acceleration  sensi- 
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ICilirfE  J7.  'Jynd;  ic  Thermal  Effect;  'J.  =  60K 
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FIGURE  38.  Dynamic  Thermal  Effect; 
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FIGURE  40.  Dynai.ic  Thermal  Effect;  ATe  =  0.5  and  1  mK 
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FIGURE  41.  Dynamic  Thermal  Effect;  aT„  =0.5 
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FIGURE  42.  Ambient  Accelerations:  Road/Rail/Ship 


tivity  for  accelerations  in  the  thickness  direction  are  the  ring-supported 
and  grooved  configurations. 20*27  Several  varieties  of  these  are  displayed 
in  Fig.  45  and  Fig.  46;  Fig.  45  gives  the  underlying  rationale  for  using 
ring-supported  resonators  when  the  acceleration  is  applied  in  the  thick¬ 
ness  direction.  The  frequency  shift  is  proportional  to  displacement  of  the 
central  membrane,  and  because  of  the  double  cantilever  edge  conditions 
this  displacement  is  smaller  than  in  conventional  simply  supported  resona¬ 
tors.  Figure  47  defines  the  parameters  pertinent  to  the  geometry  of  the 
plano-convex  ring-supported  resonator.  These  parameters  are  listed  below 
for  resonator  B-4:  (in  mm) 

$o  =  7.950;  $.j  =  6.020;  $  =  4.064; 

tQ  =  1.000;  t-  *  0.662;  t  ■  0.700; 

R  =  120. 

The  measured  frequency  constant  was  N  =  1.770  MHz  -  mm.  In  Fig.  48  is 
shown  the  experimental  frequency  shift  versus  acceleration  for  resonator 
f-4,  an  AT-cut,  plano-convex,  ring-supported  resonator  with  acceleration 
applied  along  the  X  axis.  The  frequency  shifts  are  seen  to  saturate  at 
high  g  levels,  leading  to  the  conclusion  that  the  acceleration-frequency 
sensitivity  coefficient  i;  nonlinear.  Figure  49  gives  the  corresponding 
curves  for  resonator  8-4.  The  concavo-convex  ring-supported  geometry  in 
Fig.  46  is  a  preferred  structure  for  unidirectional  acceleration  loads 
directed  toward  the  radii  of  curvature.  The  resonator  center  is  thicker 
than  the  junction  to  the  ring  to  permit  energy  trapping. 

ACOUSTIC  VISCOSITY 


One  of  the  foremost  reasons  for  using  quartz  as  a  vibrator  material  is 
its  extremely  low  acoustic  attenuation,  as  manifested  in  the  large  values 
of  quality  factor  (Q)  measured,  typically  as  much  as-  9 

Q  =  107  /  f0,  (20) 


at  room  temperature,  where  f  is  the  operating  frequency  in  MHz.  This  is 
shown  in  Fig.  50.  The  values  in  (20)  are  met,  and  exceeded,  only  by  very 
careful  design  of  the  mounting  structure  at  frequencies  below  a  few  tens 
of  MHz29;  above  these  frequencies  the  intrinsic  internal  friction  of  the 


material  predominates.  Calculation  of 
measurement  of  the  elements  of  the 

For  the  ro ta ted- Y- cuts  of  quartz, 
mode  has  an  effective  viscosity  of 

n'  =n66  co:.2e  +n44  sinze 

Substituting  the  elastic  stiffness  c .  ■ 

■  J 


the  intrinsic  loss  follows  from  the 
acoustic  viscosity  tensor.29 

the  piezoelectrical ly  driven  shear 

+  2m4  sine  cose.  (21) 

for  the  viscosities  in  (21)  yields 
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the  correspond! ng  effective  stiffness  c'  for  this  pure  mode.  Then  tne 
motional  time  constant  .  representing  the  intrinsic  losses  is 

l;  =  /c\  (22) 

Since  . .  1 s  al so  equa 1  to 

,  --  R;  C;  =  1  / (  .  Q),  (23; 

both  the  intrinsic  0  at  an  operating  frequency  f_  -  /2tt  and  the  motional 
resistance  R:  of  the  equivalent  electrical  network  may  be  computed,  as  the 
motional  capacitance  is  separately  calculated. 

Wnen  mounting  losses  cannot  be  neglected,  it  is  found  that  these  loss¬ 
es  depend  upon  the  plate  diameter-to-thickness  ratio  r  in  the  following 
manner: 


fl  +  U/t0)?)  ,  (24) 

where  .  is  the  intrinsic  value  of  x,,  and  depends  on  the  material  used, 
being  10  to  15  femtoseconds  for  room  temperature  quartz,  and  cjQ  is  a  func¬ 
tion  of  the  resonator  design,  with  =  60  a  good  average.  This  is  demon¬ 
strated  by  Fig.  51  and  Fig.  52,  where  the  data  of  Fig.  50  and  additional 
data  are  plotted  against  c .  It  is  seen  that  for  r,  '  c  the  loss,  as 
reflected  in  x,  is  governed  by  a  curve  with  -2  slope;  beyond  t  ,  the  curve 
is  constant,  with  zero  slope.  The  calculation  of  n  has  been  expended  to 
the  general  doubly  rotated  thickness  mode  case.14 

COSMIC  RAYS 

Cosmic  rays,  and  other  unshieldable  radiation  sources,  produce  cumu¬ 
lative  contributions  to  long-term  aging.  Knowing  the  radiation  sensitivity 
of  natural  and  cultured  quartz  under  a  variety  of  treatments  permits  an 
evaluation  of  resonator  performance  due  to  this  source. 

\> 

The  raw  data  are  contained  in  Fig.  53,  Fig.  54  and  Fig  55.  Using  the 
values  for  the  photon  flux  given  and  a  number  of  simplifying  assumptions, 
the  frequency  shift  due  to  this  source  is  given  in  Fig.  56  for  two  values 
of  quartz  sensitivity.  The  unit  "rad"  is  defined  as  "The  standard  unit  of 
absorbed  dose,  equal  to  energy  absorption  of  100  ergs  per  gram  (0.01  joules 
per  kilogram)."40  The  assumptions  used  are  (a)  unprotected  crystal,  (b) 
total  incident  energy  is  absorbed,  and  (c)  sea  level  flux.'Jp  is  the  photon 
flux  in  photons  /  cm;'  -  sec.  The  fractional  frequency  shift  Af/f  is  spec¬ 
ified  as  10"s;  the  time  r ( s )  for  this  shift  to  occur  is 

t(s)  =  5  x  1012_s  years.  (25) 

Similarly  for  neutrons,  assumptions  (a)  and  (b)  listed  above,  plus 
the  sensitivity  value  Af/f  ^10-21  cm2/n.  lead  to  the  order-of  magnitude 
figures  given  in  Fig.  57  for  three  values  of  altitude.  J  n  is  the  neutron 
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COISTMT.r,.  FOR  AT- CUT  QUARTZ  PLATES  AS  A  FUICTIOI  OF  DIAMETER-  THICRIESS 
OA,  AVERAGE  VALUES 


5  HAFNER  11966) 


52.  TIME  CONSTANT,  Ti,  FOR  AT-CUT  QUARTZ  PLATES  IN 
VACUUM  AS  A  FUNCTION  OF  DIAMETER-THICKNESS  RATIO,  <pafo . 


F 


FIGURE  53.  Cosmic  Ray  Flux  versus  Particle  Energy  at  Sea  Level 
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FIGURE 


Flux/cm'  -MeV  -Sec 


- 10  Km  Altitude  (300  g/cm?) 

----  3  Km  Altitude  (720  g'cmr) 


Pedicle  Energy  (eV) 


FIGURE  55.  Flux  versus  Energy  at  3  and  10  km  Height 


FIGURE  56.  Frequency  Shift  for  Photons 
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flux  in  neutron-sec/cm2.  The  fractional  frequency  shift  is  specified  as 
1CTS,  and  the  time  t(s)  to  reach  this  shift  is 

-[ ( s )  -  101 4-s/  4?  years.  (26) 

n 

CONCLUSION 

Over  the  years,  the  stabilities  of  precision  quartz  oscillators  have 
improved  according  to  a  simple  law.  With  frequency  stability  written  as 
Af/f  =  10~s,  s  is  given  by  s  =  6  +  O.l(Y-X);  Y  is  the  year  date  and  x  = 
1920  for  large-scale  production  models;  this  is  seen  in  Fig.  58.  Labora¬ 
tory  oscillators,  therefore,  can  be  projected  to  have  stabilities  in  the 
10~14  range  in  the  time  frame  2000  if  the  current  rate  of  improvement  is 
maintained.  There  do  not  appear  to  be  intrinsic  impediments  to  achieving 
10" 14  stabilities  with  quartz  elements  by  2000,  but  in  order  to  do  so  the 
various  "subtle"  effects  affecting  stability  will  have  to  be  subjected  to 
increasingly  fine  scrutiny.  Figure  59  gives  a  broad  view  of  technology 
progression  for  crystal  resonators  from  1945  onward.42  The  most  interest¬ 
ing  developments  are  yet  to  come. 


FIGURE  58.  Frequency  Stability  versus  Yea 


FIGURE  59.  Progression  in  Crystal  Technology 
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